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CENTRAL LIMIT THEOREMS FOR SEQUENCES OF MULTIPLE 
STOCHASTIC INTEGRALS 

By David Nualart and Giovanni Peccati 

Universitat de Barcelona and Universite de Paris VI 

We characterize the convergence in distribution to a standard 
normal law for a sequence of multiple stochastic integrals of a fixed 
order with variance converging to 1. Some applications are given, in 
particular to study the limiting behavior of quadratic functionals of 
Gaussian processes. 

1. Introduction. In this paper, we characterize the convergence in dis¬ 
tribution to a normal 1V(0,1) law for a sequence of random variables 
belonging to a fixed Wiener chaos, and with variance tending to 1. We show 
that a necessary and sufficient condition for this convergence is that the 
moment of fourth order of converges to 3. If is a multiple stochas¬ 
tic integral of order n of a symmetric and square integrable kernel fk, for 
instance on [0, l] n , another necessary and sufficient condition for the above 
convergence is that, for all p = 1 ,... ,n — 1 , the contractions of order p (de¬ 
fined by f® p : = fk® P fk) converge to zero in L 2 {[ 0, l] 2 ( n- f)) as k tends to 
infinity. 

In general, we call “central limit theorem” (CLT in the sequel) any re¬ 
sult describing the weak convergence of a sequence of nonlinear functionals 
of a Gaussian process (or of a Gaussian measure) toward a standard nor¬ 
mal law. The reader is referred to Major (1981), Maruyama (1982, 1985), 
Giraitis and Surgailis (1985) and the references therein for results in this 
direction. Here, we shall observe that, in the above quoted references, the 
authors establish sufficient conditions to have a CLT in the general case of 
sequences of functionals having a possibly infinite Wiener-Ito expansion. A 
common technique used in such a study is the method of moments [see, e.g., 
Maruyama (1985)], requiring a determination of all moments associated to 
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a given functional, usually estimated by means of the so-called diagram for¬ 
mulae [see Surgailis (2000) for a detailed survey]. On the other hand, our 
techniques (which are mainly based on a stochastic calculus result due to 
Darnbis, Dubins and Schwarz [see Revuz and Yor (1999), Chapter V and 
Section 3.1]) naturally bring the need to estimate and control expressions 
related uniquely to the fourth moment of each element of the sequence F 
To this end, we apply extensively some version of the product formula for 
multiple stochastic integrals, such as the one presented in Nualart [(1995), 
Proposition 1.1.3], and perform calculations that are very close in spirit to 
the ones contained in the first part of Usttinel and Zakai (1989). 

Our results are specifically motivated by recent works on limit theorems 
for quadratic functionals of Brownian motion and Brownian bridge [see 
Deheuvels and Martynov (2004) and Peccati and Yor (2004a, b)], as well 
as Brownian sheet and related processes [see Deheuvels, Peccati and Yor 
(2004)]. We provide examples and applications, mainly related to quadratic 
functionals of a fractional Brownian motion, with Hurst parameter H > 
and of a standard Brownian sheet. 

2. The main result. Consider a separable Hilbert space H. Let {e*, : k > 
1} be a complete orthonormal system in H. For every n > 1, we denote 
H &n the nth symmetric tensor product of H. For p = 0,... ,n, and for every 
/ e H &n , we define the contraction of / of order p to be the element of 
H® 2 ( n ~p) defined by 

OO 

f 0p = (f,eii®---®ei p ) H < S p®(f,ei 1 ®---®ei p ) H ® P , 

and we denote by ( f 0p ) s its symmetrization. 

In what follows, we will write 

X = {X(h):heH} 

for an isonormal Gaussian process on H. This means that X is a centered 
Gaussian family indexed by the elements of H, defined on some probability 
space P) and such that, for every h,h' £ H, 

E (X(h)X(ti)) = (MV 

For every n > 1, we will denote by I* the isometry between H 0n equipped 
with the norm \/n!|| • ||#®n and the nth Wiener chaos of X. 

In the particular case where H = L 2 (A, A , //), (A, A) is a measurable space 
and p, is a cr-finite and nonatomic measure, then H 0n = L 2 (A n ,A® n ,n® n ) is 
the space of symmetric and square integrable functions on A n and for every 
/ £ H 0n , Iif if) is the multiple Wiener-ltd integral (of order n) of / with 
respect to X, as defined, for example, in Nualart [(1995), Section 1.1.2], 

Our main result is the following. 
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Theorem 1. Let the above notation and assumptions prevail, and fix 
n > 2. Then, for any sequence of elements {f k ■ k > 1} such that fk £ H >n 
for every k, and 

(1) lim n!||/ fc ||^®„ = lim E[I*{f k ) 2 } = 1, 

k—>+oo «—►+oo 

the following conditions are equivalent: 

(i) lim fc _> +00 E[/,f (/ fc ) 4 ] =3; 

(ii) for every p = l,...,n-l, lim fc _» +0O ||/f p ||^ 02(n _ p) = 0; 

(iii) as k goes to infinity, the sequence {I* (fk) '■ k > 1} converges in dis¬ 
tribution to a standard Gaussian random variable. 

As a consequence of this result we obtain 

Corollary 2. Fix n> 2 and f £ H® n such that E [1^ (f) 2 ] = 1- Then 
the distribution of L*(f) cannot be normal and E[/;7(/) 4 ] / 3. 

Proof. If I^if) had a normal distribution or Wn(f) A ] = 3, then, ac¬ 
cording to Theorem 1, for every p = 1,... ,n — 1, we would have /® p = 0. 
Thus, for each v £ 7P x, ( n_ P) we obtain 

0 = (f JP , V ® V) H ®2 (n-p) = || (f ,V) n<S(n-p) ||#®p, 

which implies / = 0. □ 

Remark 1. The fact that a random variable with the form I* (/), n > 1, 
cannot be Gaussian is well known. See, for example, the discussion contained 
in Janson [(1997), Chapter VI]. 

To prove Theorem 1 in its most general form, we shall first deal with 
the case of X being the Gaussian family generated by a standard Brownian 
motion on [0,1]. To this end, some further notation is needed. 

2.1. The Brownian case. For every n > 1 and T > 0, given a permutation 
7T of (1,..., n), we set 

K, T = {(ti,...,t n )eR n :T> 4 (1) > • • ■ > 4 (n) > 0}, 
and we define 

t T := = {(^1 > • • • 5 ^n) £ : T > t\ > • • • > t n > 0} 

to be the simplex contained in [0,T] n . Now, for a given n > 1, take / £ 
L 2 (Alf, dti ■ ■ ■ dt n ) = L 2 (Atf). To such an / we associate the symmetric func¬ 
tion on [0,T] n , 

f (ti, ... ,t n ) ^ ] f (^7r(l) i • • • i tir(ri )) 1 A" T (^1 > • • • 1 1n): 
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where it runs over all permutations of (1,..., n). For p = 0,..., n, the con¬ 
traction of / of order p on [0, T] is the application 

(tl, . . . ,t2n-2p) ■ ■ ■ ,hn-2p) 

= / /(^l) • • • j 'U'pi ^1) • • • i tn—p) 

J[0,T]p 

X f iT'l » • • • j ^pi tn—p+1 j ■ • • i t2n—2p) du\ ■ ■ ■ dUp 

and we note f® p,T (-) s the symmetrization of /® p,T . 

In this section, we note W = {Wf,: £ [0,1]} a standard Brownian mo¬ 
tion initialized at zero. Note that, in the terminology of the previous para¬ 
graph, the centered Gaussian space generated by W can be identified with 
an isonormal Gaussian process on H = L 2 ([0,1], dt). For every n, for every 
/ 6 L 2 (A”) and every t £ (0,1], we put 

rt rs n -1 

Jn(f)= dW Sl -- dW s J( Sl ,...,S n ) 

Jo Jo 

and also /£(/) = n!J*(/), so that 

where X is once again the isonormal process generated by W. The following 
result translates the content of Theorem 1 in the context of this section. 


Proposition 3. Let the above notation prevail, fix n> 2 and consider 
a collection {gk : k > 1} of elements of L 2 ( A”) such that 

( 2 ) lim \\g k \\ L 2 {A ny.= lim ||^*||a? = 1- 

k —>+oo 1 k —>-)-oo 


Then, the following conditions are equivalent: 

(i) lim fc ^ +00 E[J^(fi( fc ) 4 ] =3; 

(ii) for every p = l,...,n-l, lim fc ^ +0O \\gf pA \\f 01] 2 n - 2 P = 0; 

(iii) as k —► oo, the sequence {J\(gk) -k > 1} converges in distribution to 
a standard Gaussian random variable; 

(iv) for every p = 0,..., n — 2, 


lim 


j A 2(n.-l- P ) ^ Sl ' " dS2( n -l-p) dtg k P t (si,...,S 2 (n—1—p)) s ^ 


( 3 ) 

where, for every fixed t £ [0,1], gk,t stands for the function on A™ 1 
by 


= 0 , 


given 


(si,••• , S 71 —1) 1 * 9k (t, Si, . . . , S n _l). 
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Proof. We will prove the following implications: 


(iv) => (hi) =► (i) => (ii) => (iv). 


(iv) =>• (iii). Suppose that (3) holds. According to the Dambis Dubins- 
Schwarz theorem [see Revuz and Yor (1999), Chapter V], for every k there 
exists a standard Brownian motion W^ such that 


( 4 ) 


r(k) 

Jo^ftrM) 2 


Jn(9k) = W\7 . .=W_ 


(k) 


[(«- 1 )!]- 2 / 0 1 dtp‘_i(3M)) 2 ' 


We can expand the stochastic time-change in (4) by means of a multi¬ 
plication formula for multiple Wiener integrals [see, e.g., Nualart (1995), 
Proposition 1.1.3] to obtain that 



where c n _i iP =p!( n j) 1 ) 2 , and therefore 

(r-s Jo dt(In-l(gk,t )) 2 _ I, m 2 Sp=0 C n-l,p/o ftt ( J 2(n-l-p) r (Os)) 
(5) [(n — l)!] 2 [(n — 1) !] 2 ' 

Moreover, thanks to a stochastic Pubini theorem for multiple integrals, 
we obtain that for any p = 0,..., n — 2, 

/o^Vi-p)^-)*)) 

[2(n— 1 — p)]! 

= dW « ‘ {Jl I' • • • ■ »*.-!-,)>.) • 

Now take the L 2 -norm of the right-hand side of (5), let k go to infinity 
and observe that if (2) and (3) are verified, then the pair 

^w!~ k \ [(n — l)!] -2 dt (^-i^fc.t)) 2 ) 

converges weakly to (W., 1), so that the conclusion is immediately achieved 
by using, for example, Theorem 3.1 in Whitt (1980), as well as formula (4). 
(iii) => (i). Trivial, given condition (2). 

(i) =>• (ii). Observe first that, for every k, 

E[4( 5fc ) 4 ] = (n!)- 4 E[^fe) 4 ]. 
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Now, as 

n / \ 2 

In(~9k) 2 = X> ! (”) lln-2 P {9t V '\-)s) 


p =0 


n -1 , s 2 

= ™'ll9*llfo,i]- + ■'(•).) + E ( ” 4-2p(9^ (•).), 

p=l Vp / 

due again to the multiplication formula, we obtain immediately 

E Kite) 1 ] = llallii + jf^ll9?"' I 0),llfo, 1] .» 


+ (» 0 ' J E 

P=1 L 


n —1 r / \ 2 - 

, ( n 


p\ 


(2n-2p)!||^’ 1 (.) s |l[o,i ] — 


Now call Il 2 n the set of the (2re)! permutations of the set (l,...,2n), 
whose generic elements will be denoted by 7r, n 1 , and so on. On such a set, 
we introduce the following notation: for p = 0,...,n, we write 


p / 
7T ~ 7T 


if the set (7r(l),..., 7r(n)) 0 (V(l),..., vr'(n)) contains exactly p elements. 
Note that, for a given 7r G Il 2 n , there are (”)"(re!) 2 permutations it' such 
that 


/ p ^ 

TT ~ 7T. 


Moreover, it is easily seen that if n' or 7T / ~7 t, then 

J 0 dai • • • da 2 n g k (^7r(l); • • • ?^7r(2n))l/fc (^ , 7T' / (1) ? • • • 5 ^'7r / (2n)) 


= \\h\\lo,i]n = (n\) 2 h k \\ 4 K , 


so that 


(2n)! 

(n!) 


4 nyfc 0 ’ (')sII[o,i] 2n 


(n!) 4 (2n)! 


E 

7ren 2 „ 


Ikllli E ( n! ) 2 + E ( n! ) 2 


. ,o 

7T ~ 7T 


n —1 


10 - X p 

+ E E / dai---da2nfff°’ 1 (a 7r )5f°’ 1 (a7T') 

•do,!] 2 ’'* 

^ 7T'~ 7T 

where = (a^),..., a^n )) ■ Since > for p = 1 ,..., n - 1 , 

E / dai-"da2„flf°’ 1 (a T )5f 0 ’ 1 (a 7r /) = f") (n!) 2 ||^® p ' 1 ||f 0)1] 

p L* 1 , 1 J \ Jr / 


2n—2p 5 









SEQUENCES OF MULTIPLE STOCHASTIC INTEGRALS 


7 


we get 


(2 n)! 


n- 1 ||n®P’ 1 ||2 

k l![ 0 ,l] 2n 2p 


), 1 m 1(2 _oil „ 114 , "[ 0 ,lj 2n -' 

(^ 7)4 ll«fc ws||[ 0 ,l] 2 « - 2 ||<?fc|| A l + 2 ^ ( p |( n _ p )|)2 

and therefore 

E[4( 5fe ) 4 ]=3|| 5fc || 4 Ai 


n—1 

+ ^(p!(n-p)!) 

p=i 


-2 


Ik II[0,1] 2 "- 2 p + 


2n — 2p 
n — p 


\ ||2 

Ik W s II [O,!] 2 "— 2p 


This yields in particular that, if E[J 4 (< 7 fc) 4 ] converges to 3, and (2) holds, 
then necessarily, for every p = 1,... ,n — 1, 


( 6 ) 


1 ™ll^f P,1 ||fo I i] 2 «-?p =°- 


(ii) => (iv). We introduce some notation: for any m > 1, x m is shorthand 
for a vector (xi,..., x m ) G ]R m , x m = maxj(xj), and dx m stands for Lebesgue 
measure on M m . To conclude the proof, we shall show that, for p = 1 ,,n — 
1, condition (6) implies necessarily that 

2 


lim 

k 


[ ds n —p f dr n - p ( f dtcjfoj. (s n —p,T n —p)\ —0. 

J[0,l] n ~ p J[0,l} n ~P S n —p\/ T n—p ’ / 


^ Sn-pVT n _p 

Now, since for every (t, s 2 , • ■ •, s n ) G [0, l] n , 

9k,t{s 2 , . . • , S n ) = l[Q^]n-l (s 2 : • • • j S n )gk(t, S2 , • • • , S n ), 

we obtain that 

[ ds n _p f d,T n —p( f dtg^P (s n _p, 

J[0,ll"-P J[0,ll"-P \Js n _ p V-r„_ p ’ / 


= f 1 dt! [ dvp_i f 1 dt [ 

Jo J [ 0 ,P]p _ 1 do d[o,t]p _1 


(iu p-l l(Up_l<t,Vp_l<t') 


( / ds n —pgk (f,,s n _p)p/j(i ,Vp_i,s n _p) ) 

V*/ [O.tAt'l ’ 1-p / 


< / dt 1 


dv p _1 / dt 


/o dfo,!]?- 1 


du T ,_ 


/o dfo,!]?- 1 


p-l 


( / ds n —p c/h (t, Up—1, s n —p)gk(t , Vp_i,s n—p) ) 

\dfO,tAPl n -P / 


X 
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'[0,1]"-p 

r-1 


dSn—p 


X 


[0,l] n_p 
dt 


(j'T 
• n—p 


O dt f d\lp—i Qk{t^ Up_x, S n _p)(t, Up_i, T n —p^\ . 

s n _ p Vf n _p J[ 0,1]p _1 / 


Now, since condition (6) holds by assumption, we know that 


m = [ 

J[o,iy 


dSn—p 


[ 0 , 1 ]" 


(]'~r 

Uj 1 n—p 


( / d\lp Pfc(llp, (tip, T n _p) j 

V[o,i]p / 


converges to zero, as k goes to infinity. We may expand the above expression 
and write 


m = f 


dSn—p 


[0,l] n -P 

Sn-pV n—p 


[ 0 , 1 ]’' 


d'T 

Uj i n—p 


X 


+ 


dt 


o 


'[o,i] n -p 


/7c 

u^n-p 


/ 

/[o,i] n -p 


dllp— 1 CJk (t, Up—1 5 S n—p)9k (^j Up_i, Tn—p) 


r]^r 
1 n—p 


O dt f p/j(t, u.p_i, s n _p)p/ c (t, u.p_i, 

Sn-pVTn—p J[ O.lJP- 1 / 


+ 2 


'[ 0 , 1 ]” 


dSn—p 


[0,1]” 


d'T 

u 1 n—p 


Sn-p V n—p 


dt 


'[o,i]p- 


dUp_l C/ k (t, Up__l, S n —p)g k (t, Up_l, T n—p) 


X f [ dt' [ dv p -ig k (t',Vp-i,s n -p)g k (t',Vp-i,T n -p)) 

\2s n _pV-r n _p J[0,1]P _1 / 

= A\{k) + A2(k) + A^(k). 

In particular, the Pubini theorem yields 

A;s(k) = / dtdup_i / dt'dvp-i 
J[o,i]p 7[o,i]p 

x ^ ds n —pl.(j;i/\t<g n _ p <ti\/t')g k (t,Up—i,s n —p')g k (t,'Vp—i,s n —p) S J 

+ 2 / dt dup-i / dt'dv v _i 
jfo.ilp 2ro,ilp 


[ 0 , 1 ]” 


dSn—p l(t'At<s„_p<t'Vi) 
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X 9k (ti Up_i, S n—p)9k{t , Vp—l, S n _p)^ 

( / d'Tn—pl(T- 7l - v <t , /\t)9k(t,Up—i,Tn—p)gk{t,' v p—i,Tn—p))’ 


Since 


A 2 (k) = / dtdup-i / dt' 

J[o,i] p J[o,i]p 

X I / dT n —p l(-r <t'/\t) 9 k (tj Up_i, T n —p)g k (t , Vp_i, T n —p) I , 

\J[ 0 ,l] n - p ” / 

we obtain that, for every k, A 2 (k) + A£k) > 0, and therefore that Ai(fc) 
must converge to zero as k tends to infinity. This gives immediately the 
desired conclusion. □ 

Remark 2. There exists an elegant explanation of the implication (iii) 
=>■ (iv) for the case n = 2. To do this, consider for simplicity a family 
{lK9k)}k>i of multiple integrals with symmetric kernels {gk}k> 1 such that 
X ]2 = 1. Then, by using the formula 


E[exp(iA/ 2 ( 5 fc ))] = exp 




Lj>2 


where is the jth power of the Hilbert-Schmidt operator associated to the 
kernel g k (■,■), it is easy to show that if (9k) converges in law to a standard 
Gaussian random variable as k goes to infinity, then necessarily 

limTr(^) = 0 

k 

for every j > 2. In particular, by taking j = 4, we obtain that 

limTr(g^) = lim [ dt [ du( [ g k (t,r)g k (u,r) dr) =0. 
k k Jo Jo \J 0 / 

But some calculations (analogous to the ones at the end of the proof of 
Proposition 3) yield 


n l r 1 / r 1 

dt du[ g k (t,r)g k (u,r)dr 
1 Jo \J 0 


rtVu 


= I dt I dui I g k (t,r)g k (u,r) dr 

Jo Jo \J 0 

r 1 rl / rr/\r' 

+ dr dr'( / dt g k (t,r)g k (t,r' 

Jo Jo \J 0 
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r 1 r 1 / rr\Zr' \ 2 

+ dr dr'l / dt g k (t,r)g k (t,r')\ 

JO Jo \JrAr' / 

/*1 /*1 / /*rVr / \ / rrAr' \ 

+ 2/ dr dr'l dt g k (t,r)g k (t,r')\l dt g k (t,r)g k (t,r')), 

Jo Jo \Jrf\r' J \J 0 / 

thus giving the desired conclusion. 

A slight refinement of Proposition 3 is the following. 

Proposition 4. Let {Jn{g k ): k > 1} be a sequence of iterated integrals 
as in the statement of Proposition 3 satisfying either condition (i), (ii), (iii) 
or (iv). Suppose, moreover, that the sequence 

dtlA^gk^-.kPl^ 

converges to zero in probability as k goes to infinity. Then, the pair 

C Jn(9k),W) 

converges in distribution to 

(N, W), 


where N is a standard Gaussian random variable independent of W. 


Proof. We shall only prove the asymptotic independence, which in this 
case is given by an application of an asymptotic version of Knight’s theorem, 
such as the one stated, for instance, in Revuz and Yor [(1999), Chapter XIII], 
□ 


2.2. Proof of Theorem 1. Let the assumptions and notation of Theorem 
1 prevail. Since H is a separable Hilbert space, for every n > 1 there exists an 
application i n (-) from H &n onto L 2 S {[ 0, l] n , dt\ ■ ■ • dt n ), such that, for every n, 
i n {') is an isometry and, moreover, the following equality holds for every 
/ € H Qn : 

n\4(i n (f)) = n\ j 1 --- [ tn ~ 1 dW tl ---dW tn i n (f)(h,...,t n f= ) lX(f ), 

Jo Jo 

where IP is a standard Brownian motion on [0,1]. It is therefore clear that 
the sequence {f k : k > 1} in the statement of Theorem 1 is such that: 

1. lim fc ^. +00 n!||z n (/ fe )||| )il ]„ = 1. 

2. lim^E [In(fk) 4 ] = 3 if, and only if, 

(n!) 4 limE[J 4 (i n (/ fc )) 4 ] = 3. 

k 
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3. For every p = 0,...,n-l, lim fc \\f® p \\ 2 He2 „- 2 p = 0 if, and only if, 

lim ||in(/fc)® p ||f 0i i]2n-2 P = 0. 

4. For k —> +oo, the convergence 


ei/j'Jsv, 


where N is a standard Gaussian random variable, takes place if, and only 

if, 

The proof is easily concluded. 


3. Examples and applications. 


3.1. Quadratic functionals of a fractional Brownian motion. Consider a 
fractional Brownian motion (fBm) B H = {B^,t G [0,1]} of Hurst parameter 
H G (0,1). That is, B H is a centered Gaussian process with the covariance 
function 

(7) R H (t,s) = l(s 2H + t 2H -\t-s\ 2H ). 

According to the so-called Jeulin's lemma [see Jeulin (1980), Lemma 1, 
page 44], with probability 1, and for every H, 

(8) J q ( 5 f) 2 ^TT = +°°- 

In what follows, we shall use the results of the previous sections to charac¬ 
terize the speed at which the two quantities 

(9) F p = j\^{B^fdt and L £ = J\b?) 2 ^ t 

diverge to infinity, respectively, when f3 tends to —H — -f, and when e tends 
to 0. In particular, to be able to apply Theorem 1, we will mainly concentrate 
on the case H > f. 

We will first find the Wiener chaos expansion of Fp. Clearly, 

~ 2j3 + 2H + 1' 

If D denotes the derivative operator [see Alos and Nualart (2003)], we have 

D s Fp = 2 f 1 t w B? dt 

J S 
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and 


D r D s Fp = ^-y(l - (« v r) 2(S+1 ). 

;k (1987)] yiel 

h{ i-(-v-) 2/J+1 ), 


2/3 + 1 

Stroock’s formula [see Stroock (1987)] yields 
(10) Fp= 1 1 


+ 


2/3 + 237 + 1 2/3 + 1 

where /2 denotes the double stochastic integral with respect to B H . 
The following result is another CLT. 


Proposition 5. If H > \, then 

(11) (2/3 + 237 + 1) f 1 t 2f5 (B?fdt 1 ^ 1, 

Jo 3l-n 1/2 

and 

(12) (2/3 + 237 +1) -1 / 2 ^(2/3 + 237 + 1) t 2 ^(B^) 2 dt — l^j ^ ^* 1/2 ChN ’ 

where N is a random variable with distribution 3V(0,1) and c# is a constant 
depending on 37. 


Proof. We denote by £ the set of step functions on [0,1]. Let H be the 
Hilbert space defined as the closure of £ with respect to the scalar product 

( 1 [o,t], 1 [o,s])n = R H(. t , s )= aH J o J o \r -u\ 2H ~ 2 dudr, 

where an = 37(237 — 1). The mapping l[ 0 t ] —► B t can be extended to an 
isometry between TL and the first chaos of B H . 

Let us compute 

E(7 2 (l - (• v -) 2/3+1 )) 2 = 2||l - (• v -) 2/m ll n®n 

= 2 a 2 H [ (1 - (s V r) 2/3+1 )( 1 - (t V u) 20+1 ) 

J[ o,i ] 4 

x \t — s\ 2H ~ 2 \r — u\ 2H ~ 2 ds dt dr du. 

Simple computations lead to 

lirn (2/3 + 237 + 1) 2 E(3+1 - (• V -) 2/3+1 )) 2 = 0, 

M-H-l/2 

because E(3 2 (l — (s V r) 2 ^ +1 )) 2 behaves as (2/3 + 237 + l)” 1 as /3 j —37 — I. 
Therefore, (11) follows. In order to show (12), set 

G/3 = (2/3 + 237 + 1) —1/2 ((2/3 + 233 + 1 )F P - 1) 

= (2/? + 2ff + 1)1/2 3 2(1 - (• V -) 2/3+1 ). 

2/3 + 1 zk v j j 
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We have 


lim 

PI—H—1/2 


E(G}) 


= lim 2 
PI—H—1/2 


(2/3 + 2H + l)a 2 H 
(2/3+ 1) 2 


/ (l-(sVr) 2/3+1 )(l-(tV«) 2m ) 

J[ o,i ] 4 


x [t — s| 


\r — u\ 


' ds dt dr du 


a H 


lim (2/3 + 2/t+l) 


2H 2 pi-H-i/2 

x / (s Vr) 2f3+1 (i Vu) 2/3+1 
•mil 4 


x \t — s| 


r — u 


' ds dt dr du 


= c 2 h . 


Taking into account Theorem 1, it suffices to show that condition (ii) holds, 
that is, 


lim 


(2/3 + 2H + l) 2 


Pl-H- 1/2 (2/3+ 1) 4 

/ (l-(-Vr) w ) 

■^[0,1] 2 

x (1 — (• V r') 2/3+1 )|r — r'\ 2H ~ 2 dr dr' 

The above norm can be written as 

a% [ (1 - (s V r) 2/3+1 )( 1 -(tV r') 2/3+1 ) 

v ro.ii 8 


= 0 . 


WigiW 


x (1 - (s' V Ii) 2/m )(l - (t' V u ') 2f}+1 ) 
x | r _/|2H-2| u _^|2H-2 

x |s — s'| 2// “ 2 |t — t'\ 2H ~ 2 dr dr' du du' ds ds' dt dt', 


and it is of the order of (2/3 + 2H + 1) 1 as /3 | —H — i, which implies the 
desired result. □ 


Note that we have also the following noncentral limit theorem. 
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Proposition 6. For any H <E (0,1), 

(13) (2/3 + 277 + 1) f 1 t 2/3 (£f f dt (Rf ) 2 . 

Jo /3T+oo 

Proof. The convergence (13) follows from (10) and the fact that 

h( 1) = (#f ) 2 - I- □ 

Remark 3. The combination of Propositions 5 and 6 generalizes results 
previously obtained for a standard Brownian motion and a standard Brow¬ 
nian bridge. In particular, analogs of formulae (11) and (13) for the case 
H = 4 are proved in Deheuvels and Martynov (2004) and Peccati and Yor 
(2004a), whereas a version of (12) for H = | is obtained in Peccati and Yor 
(2004a). 

Next, simple calculations yield that the Wiener chaos decomposition of 
L e , for any e > 0, is given by 

L e = log ^ - 7 ^ 2(1 - (e V ■ V r 2H \ 

where the second summand on the right-hand side is the double stochastic 
integral, with respect to B H , of the symmetric function 

(s,t) 1 —s- 1 — (s V t V s)~ 2H . 

This yields the following CLT. 

Proposition 7. If H>\, then 

(14) (log)) 7 ( i »- 1 °s))S*" JV - 

where N is a standard Gaussian random variable, and ku is a real constant 
depending exclusively on H. 


Proof. We keep the notation in the proof of Proposition 5. Since for 
every t > 0 the Wiener chaos decomposition of (Bf 1 ) 2 is given by 

(B ?) 2 = t 2H + I 2 (l { .v<t)), 

it is easily verified that 

rl dt f 1 ds 


E [Li] = 


= log ■ 


t 2H + 1 

1 ' " 


S 2H+1 


E[(Rf) 2 (Rf) 2 ] 


+ ■ 


1 /'* dt 


ds 


(t 2H + s 


-2H 


2H\2 


f2H+l J S 2H+1 


2 
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and consequently that 


E[L 2 £ }-E[L £ ] 2 = E[L 2 £ }-(log^j' 


behaves like log(l/e) when e tends to zero. Thanks to Theorem 1, to prove 
(14) it is now sufficient to show that 


lim[log(l/e)] 

ej.0 


-2 


X 


f (l-('VrVe) 2i? )(l — (• V r' V e) 

J [ 0 , 1] 2 

= 0 . 

Now, it is clear that 

\-2H 


~ 2H )\r - r'\ 2H ~ 2 dr dr' 


2 

h®u 


\r — r 


l\2H-2 


[ (1- (-VrVe) 2H )(1 — (• V r' V e) 

J[ o,i ] 2 

= ai / (1 - (s V r V e)~ 2H ){\ - (tVr'V e) -2 ^) 

4fn.ll 8 


dr dr' 




x (1 - (s'VtiV e)“ )(1 - (t' Vm'v e) _2H ) 

x | r _ r /|2H-2, _ n /|2H-2,_ s /|2H-2 


X |t — t‘ 


l\2H—2 


dr dr du dv! ds ds dt dt ', 


and one can show that the latter quantity is asymptotic to log(l/e) for e 
converging to zero. This concludes the proof. □ 

Remark 4. An analog of Proposition 7 for the case H = \ is proved in 
Peccati and Yor (2004b). 


3.2. Quadratic functionals of a Brownian sheet. We now extend the re¬ 
sults of Deheuvels and Martynov (2004) and Peccati and Yor (2004a) to the 
case of a Brownian sheet W on [0, l] n , that is, 

W = {W (x \,..., x n ) : (si,... , x n ) G [0, l] n } 
is a centered Gaussian process with covariance function 

n 

E[W(aq,..., x n ) W(yi,..., y n )\ = JJ (y* A xf) . 

i =1 

Note that the Gaussian space generated by W can be identified with an 
isonormal Gaussian process on L 2 ([ 0, l] n , dx\ ■ ■ ■ dx n ). 
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In particular, we are interested in the limiting behavior of the two func¬ 
tionals, defined, respectively, for vectors (3 = (f3 \,... ,/3 n ) such that A > —1 
and for e > 0, 


An = 


B e = 


'[ 0 , 1 ]’ 


dx i • • • dx n f x 2/3i ) W(xi 


, . . . , , 


,2=1 


f dx i•• • dx r 

[e,l] n (x\ ■ ■■X n )‘‘ 


■W (xi, . . . , Xnf 


when (3 converges to (—1,...,—1) and e converges to zero. We recall that, 
due again to Jeulin’s lemma, with probability 1, 


lim An = lim B e = +oo. 

/ M - 1 ,-,- 1 ) P 6->0 


Remark 5. The law of the random variable Ag, for a fixed vector f3 
and for n = 2, is studied in detail in Deheuvels, Peccati and Yor (2004). 


Now, for every (3 and for every e, we write the Wiener chaos decomposi¬ 
tions of Ag and B e which are given by 

n n ^ / n \ 

Ag = lim + 2 )- 1 + n ——-h - 0* v Vi) Wl+l ) , 

i =i + i Vi=i / 

n((*i Vy, Ve)" 1 - 1) V 

i=1 / 

where I 2 now stands for a double Wiener integral with respect to W. Then, 
calculations analogous to those performed in the previous section yield, 
thanks to Theorem 1, another CLT. 



Proposition 8. Let the above notation and assumptions prevail. 

(i) When f3 —+ (— 1,..., — 1), 

n / n 

nm +2) i/2 Ag - nm+ 2)- 1 
2=1 V 2=1 

n / n 

= um + 2)- 1 / 2 um + 2 ) a p -1 

i =1 V 2=1 

converges in distribution to y/2™N(0, 1), where IV (0,1) indicates a standard 
Gaussian random variable, and also (Il?;=i(2A + 2)Ag — 1) converges to zero 
in L 2 . 
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(ii) When e —> 0, 



converges in distribution to \/2NV(0,1). 
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